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The two-orbital Hubbard model for the pnictides is studied numerically using the real-space 
Hartree-Fock approximation on finite clusters. Upon electron doping, states with a nonuniform 
distribution of charge are stabilized. The observed patterns correspond to charge stripes oriented 
perpendicular to the direction of the spin stripes of the undoped magnetic ground state. While 
these charge striped states are robust when the undoped state has a Hubbard gap, their existence 
when the intermediate-coupling magnetic metallic state of pnictides is doped was also observed for 
particular model parameters. Results for hole doping and implications for recent experiments that 
reported electronic nematic states and spin incommensurability are also briefly discussed. 



I. INTRODUCTION 

The recent discovery of the iron-based pnictide 
superconductors^ and the perception that the pairing 
mechanism may be triggered by antiferromagnetic fluc- 
tuations, similarly as in the cuprates, has motivated a 
considerable theoretical efFort to understand the proper- 
ties of multiorbital Hubbard models. Using mean-field, 
random phase approximation, and computational tech- 
niques, several groups have addressed the nature of the 
undoped ground state of Hubbard models and the dom- 
inant pairing mechanisms upon dopingi^^~ Among the 
most robust results recently unveiled are (i) the existence 
of an intermediate Hubbard U coupling regime where the 
undoped state is simultaneously antiferromagnetic and 
metallio^ii as found in transport and neutron scattering 
experimentSfi and (ii) the presence of superconducting 
tendencies upon charge doping mainly in the Aig chan- 
nel (extended s±, with nodes) p'"—'^ with the nodal pair- 
ing states Big and B2g located close in cnergyj^iii^ 

While the superconductivity and (tt, 0) magnetism 
have received most of the attention, recent experiments 
have revealed an even more complex behavior in pnic- 
tides. For example, inelastic neutron scattering (INS) 
experiments reported evidence for spin incommensura- 
bility (spin IC) in the superconducting state of the 
hole very overdoped pnictide KFe2As2i^ This mate- 
rial does not have electron-pockets in its band struc- 
ture due to the heavy doping (50%) and, as a conse- 
quence, nesting mechanisms cannot produce the spin IC, 
whose origin then remains puzzling. Neutron studies for 
FeSco.sTco.g^'^ have also reported low-energy spin IC 
peaks near the spin resonancCfii Electron-doped pnic- 
tides such as Ba(Fe,Co)2As2 do not show spin IC at low 
energies but above 50 meV there is a splitting of the 
(tt, 0) peak. 

In parallel to these developments, scanning tun- 
neling microscopy (STM) experiments on underdoped 
Ca(Feo.97Coo.o3)2As2 reported the existence of static 
unidirectional electronic nanostructures of dimension 8 



apc-Fc (with apc-Fo the inter-iron distance) along the 
a-axisji^ This electronic nematic order state is qualita- 
tively similar to those widely discussed in other materials, 
such as in the high-temperature superconductors based 
on copper fi^ii^ Note that the lines of spins that point in 
the same direction, namely the orientation of the {tt, 0) 
spin stripes, is actually along the 6-axis, namely the ne- 
matic order and the spin-stripe order are perpendicular 
to one another. 

The spin IC and nematic order reveal new similarities 
between the Cu- and Fe-based high temperature super- 
conductors, more than expected considering the widely 
accepted perception that the former are in the strongly 
correlated regime of Hubbard U couplings while the lat- 
ter are at weak or intermediate coupling. In the Cu-oxide 
context, both spin IC and several electronic anisotropics 
in transport have been rationalized in terms of charge 
striped states, either static or dynamicJ^ Could it be that 
similar states are also of relevance in the pnictides? 

With this motivation, in the present effort the possi- 
ble existence of charge striped states in Hubbard models 
for the pnictides will be explored. Within the real-space 
Hartree-Fock approximation, it will be shown below that 
the two-orbital Hubbard model indeed displays charge 
striped states upon doping. The charge, spin, and orbital 
properties of these states will be discussed. The proper- 
ties of the striped states reported below are robust when 
the coupling U is above the critical value needed to open 
a gap in the undoped limit. Considering the complex 
nature of multiorbital systems, it is important to docu- 
ment the properties of these striped states even in gapped 
regimes that at first sight are unrelated to undoped pnic- 
tides that are known to be (bad) metals. However, the 
charge amplitude of the stripes (namely the difference 
between the largest and smallest values of the charge at 
every site in the striped state) decreases with decreas- 
ing gap, and an interesting observation is that (weak) 
stripes are still found upon doping the intermediate cou- 
pling magnetic and metallic state of the undoped limit 
for particular hopping parameters. This result is concep- 
tually different from the previous investigations of stripes 
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in doped largc-[/ Hubbard insulators. 

II. MODEL AND METHOD 

In this publication, the Hubbard model^Ti^ based on 
the dxz and dy^ Fe orbitals will be investigated. The 
use of these two orbitals is reasonable since they provide 
the largest contribution to the pnictide's Fermi SurfaceJ^ 
Studies involving models with more orbitals are certainly 
important, but they will be addressed only in future in- 
vestigations after introducing the present novel findings 
based on two orbitals. A variety of other results that have 
recently been gathered for the two-orbital model^i^i^ have 
already clearly shown that this model is in good agree- 
ment with experiments in several respects, including the 
existence of (tt, 0)-(0, tt) magnetic order, a Fermi surface 
that by construction is similar to those found in photoe- 
mission, and pairing channels that include the Aig state 
(of which the s± pairing is a special case) as well as com- 
petitors with Big and symmetries. Note also that the 
use of five orbital models would demand a much larger 
computational effort than reported here, since the diago- 
nalizations used in the iterative procedure to solve the HF 
equations scales as (2A^)''', with N the number of sites of 
the real-space cluster and 2 the number of orbitals. Thus, 
a similar calculation for 5 orbitals will require (5/2)^'-^39 
times more computer time. It is for these reasons that 
it is reasonable to report first the analysis of the case of 
two orbitals, identifying the main tendencies, and then 
in the near future increase the number of orbitals. In 
the case of five orbitals, recent ab-initio calculations sup- 
plemented by the construction of a model Hamiltonian 
led to predictions for the values of all the couplings that 
will be used in our modelJ^ However, since not all pnic- 
tides have the same couplings, and since the calculation 
of Ref. [l3 involves approximations, in our present effort, 
and in future efforts with more orbitals, an exploration 
will be carried out varying the couplings in a reasonable 
range (mainly focussing on intermediate U values, as sug- 
gested in Refs.|^0), as opposed to fixing those couplings 
to a unique set. 

As hopping amplitudes in the tight-binding portion of 
the model, fits to band calculations^ as well as Slater- 
Koster (SK) hoppings^i^ will be used. The electronic in- 
teraction is given by the standard Hubbard term, which 
has an intraorbital repulsion U, a Hund coupling Jh, an 
interorbital repulsion U' (fixed as U'~U-2Jyi here), and a 
pair- hopping term with coupling J'=J]i^ The explicit 
model does not need to be reproduced here since it is 
well known from several previous investigationsj^"— In the 
standard Hartree-Fock (HF) approximation in real space, 
the first term with the intraorbital repulsion U becomes 

5(cL,f7Ci,a,cr')(Ci,a,-<TCi,a,-<7')]- The othcr terms involve 
HF expressions too long to be reproduced here, but they 
are standard and the reader can find the full HF Hamil- 
tonian in the supplementary information The mini- 



mization of the HF energy with respect to the various 
expectation values was carried out numerically via an it- 
erative process,— mainly employing 16x16 clusters, with 
a starting configuration for the expectation values cho- 
sen at random at every site (i.e. without biasing toward 
striped states). Several clusters smaller than 16x16 were 
also explored and size effects do not seem large. 



III. RESULTS 

A. Hoppings of Refs. [3Ji 

Let us start the description of our main results using 
the Slater-Koster hoppings introduced in Refs. A 
wide range of electronic densities (n) were numerically 
explored by varying the chemical potential /x. In most 
cases, the results of the HF model energy minimization 
were sufficiently clear that they admit a simple discus- 
sion, and the emphasis below will be on those special 
cases, such as (ri)~2.3. At these densities, several start- 
ing configurations for the unknown expectation values 
that appear in the HF model were chosen via a random 
number generator and it was observed that the iterative 
process leads to nearly identical solutions or, if this was 
not the case, to local minima solutions with a higher en- 
ergy. When different types of solutions were found, of 
course only those with the lowest energy were kept and 
analyzed. 

Following the energy minimization criterion, the 
charge patterns found in some other cases were slightly 
more complicated, with segments of stripes clearly 
formed at the local level but sometimes with these seg- 
ments not properly merged together to form nearly per- 
fectly spaced stripes. Since these states still have all the 
characteristics of stripe states, it may occur that for a 
particular set of couplings some stripe configurations do 
not fit properly in the cluster sizes used and geometri- 
cal frustration effects may lead to the partial breaking 
of the stripes. So while the focus below is on the most 
clear cases, the other states found are all "stripy" in na- 
ture, thus our results seem generic of a broad range of 
couplings and densities <2i 

Figure [1] (a) contains a typical HF state that has been 
observed repeatedly in our studies, even using several dif- 
ferent starting configurations for the iterations. It is clear 
that the charge is not uniformly distributed but it forms 
vertical stripes, breaking spontaneously rotational invari- 
ance. Of course, the 7r/2-rotated state is also a degener- 
ate solution (horizontal stripes) and the convergence to 
one or the other depends on the randomly chosen initial 
state for the iterative process. Note that these stripes 
are "weaker" than the usual Cu-oxides HF stripes in the 
sense that the charge difference An (~0.15 in Fig. [T](a)) 
between the maximum rimax and minimum nmin local 
charges is not as large as in the cuprates where An is of 
order 1. This simply arises from the values of the cou- 
pling U that were investigated that are smaller for pnic- 
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tides than cuprates, in units of the bandwidth. Also it is 
interesting to note that nmin in Fig. [1] (a) is 2.27. Thus, 
the regions in the striped state that have the less charge 
still deviate from the "undoped" limit, and they are elec- 
tron doped as the rest of the striped state, while in the 
cuprates the regions between the hole stripes have den- 
sities very close to those of the undoped insulator. From 
this perspective, our striped states should be metallic, 
in agreement with the density-of-states (not shown) that 
has a nonzero weight at the chemical potential ^. 
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FIG. 1: (color online) (a) Example of charge striped state 
found in the HP approximation to the two-orbital model, us- 
ing the hoppings of Refs. Hi, (n)=2.33, f/=0.8, Jh/U^O.25, 
and a 16x16 cluster. The size of the circles is linearly related 
to the charge, with the largest circles denoting nmax=2.42 
and the smallest nmin=2.27. Here, and in the other figures, 
full (open) circles are used when the local density is larger 
(smaller) than the average, (b) Mean value of the spin in 
the state shown in (a). Note the presence of domain walls at 
the location of the charge stripes, inserted in a mainly (0, tt) 
background. 



The expectation values of the spin at each site, defined 
as Si={J2a,a,a' 4,a,<7^'^'^"^ha,cr'), that corrcspoud to the 
striped state shown in Fig. [1] (a) are given in Fig. [1] (b). 
The spin pattern is dominated by the (0,7r) configura- 
tion, with spins parallel along one axis and antiparallel 
along the other. However, at the locations of the max- 
imum electronic densities in Fig. [T] (a), the (0,7r) spin 
order breaks locally and domain walls are formed)^ with 
the spins where the local charge is maximized presenting 
an orientation nearly perpendicular to the rest. In a HF 
minimization problem, where all the expectation values 
are entangled, it is difhcult to establish if the spin state 
with walls drives the charge stripes or vice versa. But 
by analogy with the cuprates, it can be expected that 
the spin state (0,7r) tries to expel the extra charge since 
it is disruptive to such a state, and such excess charge 
is located at the walls where the spin does not maintain 
locally the (0,7r) order. Of course, all these statements 
made for vertical stripes and spin (0, tt) order are the 
same for the rotated degenerate configuration with hori- 
zontal stripes and (tt, 0) spin order. 

The HF state Figs. [T] (a,b) has other interesting prop- 
erties: (i) The orbital dyz is more populated than d^z in 
regions with the lowest local charge, i.e. where the spin 



order is locally ^(0, tt). In our opinion, this should not be 
considered as indicative of long-range orbital order but 
instead it indicates an orbital weight "redistribution" in- 
duced by the spin order that breaks rotational invariance, 
as recently discussed elsewhere»^ Where the local charge 
is the largest, on the other hand, both orbitals are pop- 
ulated approximately equally, (ii) The deviations from a 
perfect (0,7r) spin state induce a small shift away from 
(0, tt) in the spin structure factor peak position, towards 
(7r,7r). Then, the HF striped states can produce spin 
incommensurability, as found in some experiments (see 
Introduction), (in) Note that the charge stripes and the 
spin stripes (i.e. the lines of parallel spins in the (0,7r) 
arrangement) arc perpendicular to one another, as ob- 
served in the STM experiments where nematic order was 
reported 
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FIG. 2: (color online) (a) HF charge striped state using the 
hopping amplitudes of Ilefs. [^Q, at (n)=2.45, [7=0.8, and 
JH/t/=0.25. The size of the circles is proportional to the 
charge, with nmax=2.48 and nniin=2.41. (b) Same as (a) but 
for {n)=1.83, C/=1.0, with nmax=1.84 and nmin=1.81. (c) 
(tt, 0) antiferromagnetic order parameter m, charge gap A 
(from (n) vs. /i), and An=nmax-'^min at (n)~2.3, as a function 
of U . (d) (n) vs. fj, at C/=1.0, Jyi/U=Q.2^, and various lattice 
sizes, suggesting that size effects are small in this quantity. 
Results at others C/'s appear equally well converged. 



Figure [2] (a) contains another typical HF state with 
charge stripes that was found in our studies, at an elec- 
tronic density (n) larger than in Fig. [T] (a), thus con- 
comitantly inducing a reduction in the distance between 
stripes. It is interesting to observe that An is smaller to 
that in Fig.[l](a), suggesting that there must be an "op- 
timal" doping for vertical/horizontal stripe formation, as 
found in cuprate's investigations, since in the undoped 
case (n)=2.0 there are no stripes and thus An=0. Also, 
note that our results even farther from (n)=2.0 indicate 
a variety of complex patterns, with checkerboards, diago- 
nal stripes, and other arrangements. However, since their 
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associated spin orders are very different from the (tt, 0) 
(or (0, tt)) state prevailing in pnictides, and the doping is 
too large to be compared with available experimental re- 
sults, those exotic states will not be described here. Pre- 
sumably when the electronic density deviates substan- 
tially from the undoped case, then our results reach a 
regime not yet found experimentally in these materials, 
or the two-orbital model breaks down. 

Figure [2] (b) contains results for the case of hole 
doping, and here the stripes are found to be weaker, 
with An^O.03. At least for the hopping amplitudes of 
Refs. HH, the stripe tendencies appear stronger for elec- 
tron doping than for hole doping. Also, in the hole-doped 
case the lines of parallel spins tend to run parallel to the 
charge stripes, instead of perpendicular as in the electron- 
doped case. This parallel vs. perpendicular relative pat- 
terns of charge and spin may depend on details and may 
not be universal, thus they need to be investigated in 
more details in the future. 

Figured] (c) displays An vs. U, at {n)^2.3, compared 
with results for the undoped limit, namely the (tt, 0) (or 
(0, tt)) magnetic order parameter m and the charge gap 
A. While m rapidly grows with U at Ud, A remains zero 
in a narrow region beyond Ud and only after reaching a 
second critical value Uc2 it slowly increases with increas- 
ing U. This intermediate region is both metallic and 
magnetic, and has been emphasized before as the regime 
of interest for pnictides i^ii An follows A with decreasing 
U from strong coupling but it appears to survive, albeit 
with a small value, in the intermediate magnetic /metallic 
phase. If this result would survive the introduction of 
fluctuations beyond the HF approximation, the undoped 
metallic magnetic state would admit weak charge stripes 
upon doping, a conceptually novel result since previous 
stripe efforts have focused on their presence only when 
doping a Hubbard insulator. For those insulators the ra- 
tionale for the presence of stripes was dominated by the 
fact that extra charges tend to disrupt the regular spin 
order of the undoped case. Thus, to reduce this "dam- 
age" to spin order, the extra charge is arranged forming 
stripes. Such a way to understand the stripes of the 
cuprates can be extended to pnictides even in the mag- 
netic and metallic regime: the extra charge may need 
to be accumulated in patterns, as opposed to regularly 
distributed, to minimize the energy damage to the spin 
order. Note also that size effects do not seem large here 
since some quantities are already converged on the 16x 16 
cluster (sec Fig.[2](d)). 

B. Hoppings of Ref. 

The HF results obtained using another set of hopping 
amplitudes (Ref. H) are presented in Fig. [3] In this case, 
the analysis was slightly more involved, since the conver- 
gence of the HF energy minimization was slower, with 
the iterative process for convergence sometimes ending 
in metastable states with charge inhomogeneous irregu- 
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FIG. 3: (color online) (a) Charge striped state solution of the 
HF two-orbital Hubbard model, using the hoppings of Ref. 0, 
(n)=2.36, f/=0.50, and Jh/(7=0.25. The size of the circles is 
proportional to the charge, with nniax=2.365 and nmiii=2.361 
(a very weak charge stripe), (b) Spin state associated with 
(a), (c) Same as (a) but for [7=0.90, (n)=2.14, nmax=2.141, 
and nmin=2.139 (d) Spin state associated with (c). 



lar patterns (but still displaying stripes locally). How- 
ever, in many cases the convergence to a nearly perfect 
stripe pattern was clear. Panels Figs.[3Ka) and (c) display 
the charge order in the HF states at two different (n)'s 
and two U's with a magnetic undoped state (the hopping 
unit ti of Ref. |^ is taken as 0.2 eV here, as in Ref. Q). 
In both cases, charge stripes can be identified, with An 
growing with U, as in Fig. [DJc). At [7=0.50 the charge 
gap is negligible, so these (weak) charge stripes appear 
in the metallic magnetic state. The associated spin pat- 
terns also present a striped arrangement, and the spin 
and charge stripes are perpendicular. There are some dif- 
ferences between the stripes of Figs [1] and [3l such as the 
stripe periodicity. These differences are to be expected 
since both sets of hopping parameters were constructed 
to produce similar undoped Fermi surfaces but their full 
band structures are different. However, it must be em- 
phasized that with both sets of hoppings, charge stripes 
are found in several cases and, thus, their existence is a 
qualitative conclusion of our HF investigations. 



C. Preliminary studies of the influence of Co 
doping on a striped state 

The striped states found here could be of relevance 
to explain the static nematic state found via STM>i^ To 
test this hypothesis, a randomly-chosen set of lattice sites 
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FIG. 4: (color online) (a) Ground state of the HF Hamil- 
tonian at (n)=2.56, (7=1.0, and JH/i7=0.25, using the SK 
hoppings.— In the 16 sites shown in grey, there is an on-site 
energy e=-0.85, that simulates the presence of quenched dis- 
order, such as caused by Co doping, (b) Monte Carlo results 
(equilibrated snapshot) obtained using the two-orbital spin- 
fermion model of Ref.lli with SK hoppings,^^!^ (n)=2.55, K=- 
1, Jnn=Jnnn=0.05, magnitude of the classical spins 5=1, 
temperature 0.005 eV, and 45,000 steps. The small deviation 
from a perfect stripe is caused by temperature. 



was selected to simulate the presence of Co-dopants via 
an on-site energy e, and then the minimization of the 
HF model was again carried out in that background. A 
typical result is shown in Fig. IDJa). It is observed that 
this quenched disorder effectively cuts the long stripes 
into shorter segments, and the overall state qualitatively 
resembles those found with STM. Of course, only a tun- 
ing of couplings and electronic densities may render this 
agreement quantitative (for instance producing features 
of size 8 lattice spacings, as in experiments) and that ef- 
fort is postponed for future work. Here, simply note that 
the charge and spin stripes found in our study tend to 
be perpendicular to one another, in nice agreement with 
the STM investigations. 



filling has been here discussed. Charge stripes appear at 
several electronic densities and in a broad range of cou- 
plings. The associated spin and orbital properties of the 
HF states have been discussed as well. However, the read- 
ers should be warned that the relevance of these states 
to the real pnictides is a matter that still requires further 
work. On one hand, the charge stripes found here ap- 
pear to have an associated spin incomensurability. Since 
neutron scattering results have actually unveiled a sim- 
ilar spin incommensurability in cases where nesting can 
be excluded, then it is important to search for alterna- 
tive explanations such as those presented here. In addi- 
tion, the stripes discussed in this effort can also provide 
a qualitative rationalization of the nematic order found 
in STM experiments. On the other hand, the HF striped 
states are the more stable at couplings \J where the un- 
dopcd state is an insulator, albeit with a small gap. How- 
ever, for the hopping sets investigated here, novel (weak) 
stripes were observed in the magnetic/metallic region, 
thus avoiding this conceptual problem. But this inter- 
esting observation must be confirmed with more refined 
calculations and using models with more orbitals. 

In addition, it may occur that the surface of pnictides 
presents an effective U larger than in the bulk, due to 
the reduction of the carrier's bandwidth, complicating 
the analysis. Another issue to consider is that the use 
of models with more than two orbitals may lead to an 
effective electronic population in the d^z-dyz sector that 
is larger than 2.0, and charge stripes, magnetism, and 
metallicity may coexist in a single state at couplings 
U larger than those estimated to be realistic from the 
(ri)=2.0 two-orbital model analysis. All these scenarios 
are speculative right now, and only further work can clar- 
ify the relevance of the striped states presented here to 
the physics of the pnictides. 



D. Results for the spin fermion model 

Striped states were also studied via the two-orbital 
spin-fermion model recently proposedi^ Using Monte 
Carlo techniqueS)^^ metastabilities limited the present ef- 
fort to 8 X 8 clusters. In cases where good convergence was 
achieved, stripes were also found (Fig. |4l^b)). Thus, via 
the use of two rather different models and techniques the 
presence of charge non-uniformed states was confirmed. 

IV. CONCLUSIONS 

In summary, the numerical solution of the two-orbital 
Hubbard model in the HF approximation away from half- 
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